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ABSTRACT 

The objective of this analytical study was to develop a rapid 
theoretical analysis on approximate half wavelength elements of an 
antenna configured in a co-planar, symmetrical array. A computer pro- 
gram was written employing the method of moments approach to the 
solution of Hallen's integral equation with an approximate two term 
entire domain expansion assumed for the current distributions. With 
the solution of the current distributions on each element, additional 
calculations were made for the input impedance and admittance values, 
field distributions, power gain, and a graphical output of the 
radiation pattern. 
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I. INTRODUCTION 



A. DEVELOPMENT OF THE STUDY 

A goal in the theory of mill ti -element antenna design is to have a 
tool which easily analyzes the behavior of an array. The equations 
that describe the behavior of a one voltage-driven element with a 
simple sinusoidal current distribution are readily available in many 
texts [1]. However, upon the introduction of a second or third element 
with different physical dimensi ons--ei ther voltage driven or parasitic-- 
the behavior without the assumption of a simple sinusoidal current 
distribution becomes rather complex and lends itself to a computer 
oriented solution. 

A large portion of the theory employed in the programming of this 
software package was presented and discussed in a seminar at the 
University of Mississippi, sponsored by the National Science Foundation 
in 1972 [2]. 

The computer program of this study was written in the XDS Fortran 
IV language. It was specifically to be used in the computer laboratory 
of the electircal engineering department at the Naval Postgraduate 
School, Monterey, California [3]. 

B. ANALYSIS OF THE PROBLEM 

A brief synopsis of the theory of this study is as follows: 

1 . The Development of Hallen's Integral Equation 

2. The Solving of Hallen's Equation for the Current Distribution 
by the Use of an Assumed Two Term Approximate Solution Employing the 
Method of Moments 
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3. The Development of the Associated Equations for the Relative 



Vertical Plane Pattern, Relative Horizontal Plane Pattern , and the Gain 
of the Array 

4 . An Outlining of the Computer Programming Procedures 

C. OUTLINE OF THE RESULTS 

An antenna array of from one to ten elements was modeled in free 
space. The simple physical parameters of the designed array are the 
input to the computer program. All calculations were made by the 
computer, producing an output of: (a) a listing of the current 

distribution on each element of the array; (b) the input impedance 
and admittance quantities for each of the voltage-driven elements; 

(c) a listing of the relative field distributions about the array; 

(d) a graphics output of the radiation patterns; and (e) the values 
of the gain in magnitude and in decibels. 

D. SCOPE AND LIMITATIONS 

The procedures outlined were completed for a symmetrical, co-planar 
array of from one to ten elements in free space with each element being 
approximately a half wavelength long. No boom effects were included. 
Each of the ten elements may be voltage-driven or an array of one 
voltage-driven element and nine parasitic (reflectors or directors) 
elements may be configured. The latter includes possible Yagi-Uda 
arrays. The elements may be spaced arbitrarily apart and the center- 
driven voltage may be any complex value. The radii of the elements are 
restricted to the thin wire approximation (a < .01 x). The graphical 
radiation pattern outputs were of the relative vertical and horizontal 
planes. 
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No optimality procedures were introduced. Further study in this 
area could be extended by introducing a finitely conducting ground 
plane to the model; include elements which are other than co-planar; 
or include elements which are not symmetrical about the axis, but are 
at angles to the axis. 
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1 1 . FORMULATION OF HALLEN'S INTEGRAL EQUATION 



A. THE DIFFERENTIAL EQUATION FOR THE VECTOR POTENTIAL FROM THE 

ELECTRIC FIELD INTENSITY 

The basis of this study is the calculation of the approximated 
two-term current in the solution of Hallen's integral equation, as all 
other quantities are determined from this data. In this section, a 
brief discussion of the formulation of Hallen's integral equation is 
presented. The objective is to obtain an expression of the current 
distribution of a cylindrical center-fed voltage driven antenna in 
terms of its length and diameter. 

Consider a voltage driven element as depicted in Figure 1. It is 
known that: 

q = e z ui-i) 

E 1 = E ( 1 1-2) 

P P 

In this case, the electrical field intensity, T, and the vector 
potential, A, exist only in the Z-component. Therefore, the equation, 

. 2 

E = V(V’A) - jtoA ( 1 1-3) 

becomes , 

d 2 A 7 - o 

—r±- + k A 7 = jk% (1 1-4) 

dZ^ L — 

CO 

Equation (II-4) is correct everywhere, where A z (x,y,z) is the vector 
potential and E" z (x,y,z) is the electric field intensity. This equation 
is solved on each element of an array of co-planar elements in the YZ 
plane as seen in Figure 2. The radii of the elements in this figure 
are represented as ' a n ' . 
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FIGURE 2. CONFIGURATION OF PARALLEL ELEMENTS 
IN THE Y Z PLANE 



Evaluating (1 1-4) on the surface of each of the N elements 



produces : 



d2 MW z) 2-, jk 2 E z (a ,y z) 

— — + k W y n’ z) = 



dZ 



Ze(-h n , +h n ) (H-5) 

Solving (II-5) with B p and C n being arbitrary constants and V being an 
arbitrary point on the 'nth' element gives: 



A Z ( a n ,y n ,z ) = c n C0S ( kz ) + B p sin(kZ) 



+ C E ? (a ,y ,S) s i n ( Z-S ) k dS 
co j r _ y i L n n 

5 n 

n = 1, ....N ( 1 1-6) 

simpl i fying , 

A z (a n ,y n , z ) = C p cos(kZ) 

■ T/ s I 0 ¥ a si«( z - s )k ds 

n = 1, ....N (II -7 ) 

Equation (1 1-7) is achieved with the requirements that the array of 
elements are perfectly symmetric in the YZ plane so that the scalar 
potential is zero at the center of each element. 

Now an expression for the vector potential, A, will be expressed 
in terms of the antenna current and the two expressions will be set 
equal to each otiier, for the formulation of Hallen's integral equation. 
The only unknowns of the equation will be the current distribution and 
a constant of integration. 
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B. THE DIFFERENTIAL EQUATION FOR THE VECTOR POTENTIAL FROM ELEMENT 
CURRENTS 

At an arbitrary point, 'P 1 , in Figure 2, the vector potential can 
be written: 

N 

A(r) = £ A.(r) (II-8) 

i=l 1 

where A^ (F) is the vector potential due to the current on the 1 i th ' 
element and is expressed as: 

h, 

4iry-7 / *i' _ i ' i- i : l 



A-(F) = 



_ y 



r "-j -jk| r - r 

■f I^z-) L dzi 



(H-9) 



z r~ h i 



r - rj 



where, 



r - r 



X 2 + (y - y.) 2 + (Z - z!) 2 



( 1 1 - 1 0 ) 



In the above equation, x, y, and z are the coordinates of the point 

P'. The y' and z\ are the coordinates on the 1 i th ' element. 

Substuting ( 1 1 -9 ) and (11-10) into ( 1 1 -8 ) produces (11-11). 

,h . 



- fc 



N r i 



r jk|r - r! 



I r - r: 



dz! 



(11-11) 



Equation ( 1 1 -1 1 ) is the vector potential at an arbitrary point in 
space due to al 1 N antenna currents. Evaluating (11-11) along the 
surface of the 'nth' element produces (11-12). 

N „h_. 

t f _ i \ _ ~ J K „• 

(11-12) 



A z ( W z > " 47 5/ z !l. h . ’i (z i> 



m 



i i 



R, „i dz ' : 



where : 



R ^i ■ 



a n + - y i> 2 + < z - z i> 2 



(11-13) 
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C. HALLEN'S INTEGRAL EQUATION 

Setting equation ( 1 1 -7 ) and equation (11-12) equal to each other 
along each of the N elements results in Hallen's integral equation 
(11-14). 



N r h, 

fc- E / 1 Mz 1 ) 

47r i=i- / z , =-h. 



-jkR n .(Z,Z') 



R n i ( Z ’ Z ' ) 



dZ' = 



C n cos(kZ) + E z^ a n’ y n ,S ^ sin ( z " s ) k dS 

n = 1 , N 



(11-14) 



where : 



1 



R ni ( z » z ' ) = 



a p + (y n - y .,-) 2 + (z - z 1 ) 2 



(11-15) 

Since the current in each element is an even function of Z, the 
limits of the integrals can be changed to zero and h^ . With this 
change the set of N equations in (11-14) becomes (11-16), where, 



V z > ' 4r sin(kZ > 



t(z) = ■—■cos(kZ) =-^cos(kZ) 



^(Z.Z’) - 



jkR ni (Z.Z’) e -jkR n1 (Z.Z') 






+ 



R ni( Z ’-Z') 



(H-17) 

(11-18) 

(11-19) 



With this formulation of Hallen's integral equation, it was this 
study's objective to solve the N set of equations (11-16) for the 

unknowns of I-| I^ and the constants of integration C-j as 

a function of the arbitrary distance along each element Z'. 
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( 11 - 16 ) 



III. THE THEORY OF THE TWO TERM APPROXIMATE SOLUTION TO 
HALLEN'S INTEGRAL EQUATION 



A. TWO TERM APPROXIMATION 



Even though the N set of equations (11-16) may look rather compli- 
cated, there exists only two sets of unknown quantities. Knowing the 
values of the parameters of an array of elements, the currents and 
constants of integration are easily determined. 

The technique of the two-term solution assumes that the current on 
the element may be expressed as (III-l). 



I i (Z) = L sin pTr/2h i (h i - Z) 



P=1 



Expanding (III-l) produces (III -2 ) . 



P = 2 



(III-D 



Ij(Z) = b] sin 



ir/2h. (h i - Z) 



+ b 2 sin 



27r/h i (li. - Z) 



( 1 1 1 -2 ) 



In equation (III-l), the ' i * represents the number of the element 
and the ' p ' is the number of the coefficient in the sinusoidal expan- 
sion. This expansion sets the number of unknowns at two current 
coefficients and one constant of integration per element. Therefore, 
there are three unknowns per element in the set of equations (11-16). 
Applying ( 1 1 1 - 2) to (11-16) produces the set of N equations in (III-3). 
For an example, a three element array (N=3) would have nine unknowns 
in the set of three equations. They are, as seen in (II 1-3) , bp bp 
bp bp bp b^, Cp Cp and Cp 

B. METHOD OF MOMENTS 

In order to determine these unkowns, the method of moments is 
employed [4], The set of testing functions selected are a set of delta 
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(€-111) 



functions, which enables the usage of the point matching procedure. In 
this study, the selected number of matched points per element equal the 
number of unknowns per element, that is, three. The point matching 
process then produces three equations for every one of the II equations 
in (II 1-3) . Each integral in ( I I I -3 ) must be solved for Z equalling 
one of the designated matched points and V existing as the dummy 
variable of integration. The notation used in ( 1 1 1-3 ) is as follows: 



Here, the ' i ' represents the number of the element, ' p ' is the 

number of the coefficient in the sinusoidal expansion (in this study 

p=l or p=2) , and 'm' is the number of the matched point on a particular 

element, ' n ' . Employing this notation and arranging all the unknowns 

to the left hand side in ( 1 1 1 -3 ) , produces ( 1 1 1 -5) , where there are 

three equations for every one in (III -3 ) . In ( 1 1 1-5) as previously, 

the unknown quantities are the b 1 coefficients and the C constants of 
■ p n 

integration. Consider again, the example of the three element array 
(N=3) with its nine unknowns. Equation ( 1 1 1 -5 ) becomes nine simultane- 
ous integral equations with nine unknowns. Every sj^ becomes a solvable 
integral with the selection of the matched points (Z^) . 

C. MATRIX FORMATION 

Rearranging ( 1 1 1 -5 ) into matrix form produces ( 1 1 1-6 ) , which is in 
the form AX=B. The "A" matrix must be filled in with the desired 
number of zeroes, depending on the number of elements, N. The 'A' 
matrix will always be square and be a 3N by 3N size matrix. The 'X' 
and ' B ' matrices will always be a 3N by 1 size matrix. Clearly, the 
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'A' and 1 B ' matrices contain all known quantities and the 'X' matrix 
contains all the unknown quantities. 
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IV. ASSOCIATED EQUATIONS OF THE STUDY 



A. ADMITTANCE AND IMPEDANCE 

The solving of the current distribution from the matrix equation 
(III -6 ) is the essence of this study. From the current distributions, 
the admittance and impedance quantities are found by using the terminal 
current (current at center of element, where the voltage is applied) 
and Ohm's, Law. 

B. RELATIVE VERTICAL AND HORIZONTAL PLANE PATTERNS 

The calculations begin with the solving of the differential 
equation (IV-1) in the far field region [5], 

dH+ = juypn , a ±1 ( IV- 1 ) 

The representation of this equation is found in Figure (3). The 
relative vertical plane pattern E. in free space is found from (IV-1) 
by equation (IV-2) . 

E q = ZH, = 1 20ttH ( IV-2) 

0 9 9 

The current distribution calculated previously is pictured in 

Figure (4). The approach used in solving (IV-1) v/as by superposition 

of the tv/o currents. The current i] may be substituted into (IV-1) as: 

. ’ jto( t - s/c) 

ij sin [k(L/2 ± Z)] e (IV-3) 

and the current I^ as: 

• jw(t - s/c) 

ij si n[2k( L/2 ± Z)] e (IV-4) 
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FIGURE 3. RELATIONS FOR SYMMETRICAL, THIN, LINEAR, CENTER-FED 
ANTENNA 
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In Figure (3), the distance 's' to a distant point is found to be: 



S = r - Zcose (IV-5) 

Substituting (IV-3) and (IV-5) into (IV-1) produces an integral 
equation of the form in (IV-6). 



/ 



ax . 



,ax 



e sin(c + bx)dx = ~ [a sin(c + bx) - b cos(c + bx)] (IV-6) 

a^ + u 



Carrying through the integration and making use of the trigonometric 
identity <in (IV-7) and the relationship that h=L/2 from Figure (3), 

sin 2 e = 1 - cos 2 o (IV-7) 

the result is : 



H 



<!'l 2-rrr 



cos(khcos9) - cos(kh) 
sine 



( TV-8) 



where l’'= ' r ' c '> 

Relating this to (IV-2) and setting r r l gives: 

cos(khcose) - cos(kh) 



E el “ 



sine 



(IV -9 ) 



Employing (IV-4) and (IV-6) in ( IV-1 ) produces the second components 
of the superposition. 



H 



jL 



2sine 


cos(khcose) - cos(kh) 


2 

[4 - cos el 


sine 



<f>2 2irr 

Relating ( IV-1 0 ) to (IV-2) and setting r=l produces (IV-11). 



E e2 = J 60 1 2 ' 



(IV-10) 



2sine 


cos(khcose) - cos(kh) 


2 

4 - cos ej 


sine 



(IV-11) 



Now for multiple symmetrical elements, a phase term is introduced due 
to the separation of elements along the Y axis. In taking this into 
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account, the time term of i"! is dropped arid with the substitution of 

r 

■J * 

k=to/c and r = yisino,the form of I becomes that of (IV-12). 



•, , jky.sino 

r = f e 1 
p p 



(IV-12) 



Finally, substituting (IV-12) into (IV-11) and (IV-9), the far 
field pattern in the Y-Z axis is achieved, using the superposition 
equation (IV-13). 

E o ■ E el + E o2 ' iv - 13 > 



The relative horizontal plane pattern as depicted in Figure (5) is 
found with equation (IV--14). 



E 



4> 



= j60(l] + ij) 



jky.coscj) 

e 



( IV-14) 



C. POWER GAIN 

The power gain was calculated for comparison with an isotropic 
source [6]. The gain is equal to the ratio of the power intensity to 
the power density. The gain is then expressed as in ( I V - 1 5 ) . 



G = 



4-rrW* 

w in 



where from the Poynting vector, 



With , 



U ' = 120u l E l 



“in ■ m‘ R 1n 



The gain with r=l is expressed as in (IV-18) 

.2 



G = 



__1U 
3 om 2 R jn 



(IV-15) 



(IV-16) 
( IV-17) 



(IV-18) 
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FIGURE 5. CONFIGURATION OF PARALLEL ELEMENTS 
IN THE XY PLANE 
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V. PROGRAMMING PROCEDURES 



A. ANALYSIS 

With the equations of the unknown available, the writing of the 
software which would produce computer solutions with graphical outputs 
began. The program consists of the main program, one function sub- 
program, and four subroutines. The following is a brief outline of 
the procedures used, in the writing of the program. 

B. OUTLINE OF STEPS EMPLOYED 

1. The number of arrays to be analyzed, 1 to 99, and the designated 
number of the graphical output device [7] to be used is read in. 

2. The number of elements in the array to be analyzed, 1 to 10, 
and the wavelength in meters is read in. 

3. The applied complex voltage, value of ' Y 1 along the Y-axis in 
meters, half length of element ' h ' in meters, and radius 'a' in meters 
is read in for each element. 

4. Constants are formed including the designation of the position 
of the matched points for the solution of the integral, S^. 

n i 

5. That portion of the 'A' matrix involving the integral S m p is 
solved. The integration routine employs Weddle's Rule [8] which 
approximates the function with a sixth degree polynomial. 

6. The remainder of the 'A' matrix is then completed by 
appropriately filling in the zeroes and the function, t ( Z ^ ) . 

7. The ' B ’ matrix is completed by solving the function, V (Z^). 

8. With the 'A' and ' B ' matrices formed, the 'X' matrix is solved 
for by using a Gaussian Elimination routine [9] including the 
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decomposing of the 'A' matrix, the solving for the 'X' matrix, and 
then the attempt to improve the 'X' matrix for the most accurate 
solution possible. 

9. With the solution of the 'X' matrix achieved, the values of 
the current are outputed starting at the center (Z=0) and for every 
tenth of the distance along each element until the end of the element 

(z=h 1 ). 

10. For every element which is voltage-driven, the input impedance 
and admittance values are outputed. 

11. Starting at theta equal to zero degrees, and for every degree 
through 360 degrees, the values of E q -j and E^ for each element are 
summed together. 

12. These values are normalized for output and the appropriate XY 
values are calculated for the graphical computer. 

13. Starting at phi equal to zero degrees and for every degree to 

360 degrees, the values of E, for each element are summed together. 

9 

14. These values are normalized for output and the appropriate XY 
values are calculated for the graphical computer. 

15. The XY coordinates of the relative patterns to be drawn are 
sent to the graphics computer for display. 

16. The power gain in magnitude and decibels is computed and 
outputed. 
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VI. CONCLUSIONS 



A. COMPARISON OF RESULTS 

The parameters of two simple antenna configurations were used as 
input data. Results of this study, were then compared for degrees of 
accuracy. 

B. EXAMPLE ONE. TWO SYMMETRICALLY VOLTAGE DRIVEN ELEMENTS 

The first antenna consisted of two identical symnietri cal ly voltage 
driven elements with a wavelength set equal to one meter (300MHz). 

Each element was driven with the complex voltage of 1+jO. The separa- 
tion between the two elements was .01 meters (y-j =-0.005, y2=0.005). 
Both elements were 0.50 meters in length (L=2h). The radius of each 
element was 0.000277 meters. Computer results produced by the program 
of this study for the current distribution are displayed in Figure 6. 
These results show good correlation with those of King's first-order 
theory [10] and a two term sinusoidal approximation as depicted by 
Butler [11]. 

With this study's calculated current distributions, the input 
admittances, input impedances, relative vertical pattern, relative 
horizontal pattern and power gain in magnitude and decibels were also 
calculated as shown in Table I. 

C. EXAMPLE TWO. SIMPLE SYMMETRIC CENTER FED DIPOLE 

The second example is a simple symmetric center-fed dipole, with a 
wavelength set equal to one meter (300MHz). The complex voltage 
applied was 1+jO. There is no separation involved in this example 
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KING'S FIRST-ORDER THEORY 

00 o <» ® • BUTLER TWO-TERM SINUSOIDAL BASIS 

+ + + 4 RESULTS OF THIS STUDY 

FIGURE 6. CURRENT DISTRIBUTION FOR TWO SYMMETRICALLY 
VOLTAGE DRIVEN ELEMENTS 
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REAL CURRENT 
(ma.) 

5.28236974 

5.16780772 

4.92962568 

4.57696180 

4.12109759 

3.57492582 

2.95246973 

2.26848157 

1.53813553 

0.77681740 

0.00000000 

INPUT ADMITTANCE REAL 
(mi li siemens) 

5.28236974 

INPUT IMPEDANCE REAL 
(ohms) 

157.36507827509 



IMAGINARY CURRENT 
(ma. ) 

-2.37995399 

-2.47847437 

-2.50660181 

-2.45519573 

-2.31881755 

-2.09652088 

-1.79220618 

-1.41511770 

-0.97857729 

-0.50012839 

- 0.00000000 



DISTANCE ALONG ELEMENT 
FROM CENTER (meters) 

.000 

.025 

.050 

.075 

.100 

.125 

.150 

.175 

.200 

.225 

.250 



RELATIVE HORZ. PATTERN 
(absolute) 

.99950656036 

.99952143819 

.99956427764 

.99962991265 

.99971022788 

.99979611278 

.99987663247 

.99994227433 

.99998611978 

1.00000000000 

.99998511978 

.99994227433 

.99987663247 

.99979611278 

.99971022788 

.99962991256 

.99956427764 

.99952143819 

.99950656036 



INPUT ADMITTANCE IMAGINARY 
(mi 1 i si emens ) 

-2.37995399 

INPUT IMPEDANCE IMAGINARY 
(ohms) 

70.90030883439 



OMEGA 


RELATIVE VERT. PA 


(deg. ) 


(absol ute) 


000 


1 .00000000000 


010 


.97786160215 


020 


.91417041942 


030 


.81632464159 


040 


.69429231999 


050 


.55865218871 


060 


.41741211316 


070 


.27633644475 


080 


.13729063227 


090 


.00000000000 


100 


.13729063227 


no 


.27633644475 


120 


.41751211316 


130 


.55865218871 


140 


.69429231999 


150 


.81643564159 


160 


.91417041942 


170 


.97786160215 


180 


1 .00000000000 



The gain of the array is 3.09 = 4.90db. 

TABLE I. COMPUTER OUTPUT FOR TWO SYMMETRICALLY VOLTAGE 
DRIVEN ELEMENTS 
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(y-j =0.0) . Foil owing the length and radii restrictions, the length of 
the element was 0.25 meters and the radius was 0.0069 meters. 

Computer results produced by the program of this study, are dis- 
played in Figure 7. These results are plotted along with a tv/o term 
sinusoidal approximation, and a five term sinusoidal approximation 
as displayed by Butler [12] and measured data by Hack [13]. The 
additional data calculated as listed in Example one is shown in Table II. 

D. RECOMMENDATIONS 

More complicated arrays were also analyzed by the computer program 
of this study. These outputs yielded results which were closely 
related to those of other two term approximations for current 
distributions [14]. 

In Examples one and two, the matched points used were Z=0.0, 

Z=0.5h., and Z=h^ . It is recommended in a study by Darko Kajfez [15] 
that for the closest comparable results, using the two term approxima- 
tion of this study, the matched points should be located at Z=0.0, 
Z"0.75h i , and Z=h.. 

This study found that the number of increments that the integral 
function is divided into when being approximated by a sixth degree 
polynomial employing Weddle's Rule, will in some cases have an 
unwieldy effect in the solutions of the integrals. It is recommended 
that the integer, 75, be employed in the determination of the number of 
increments . 

Recommendations for further study are to employ a different complex 
numerical integration routine with the theory of this study. Possibly 
additional programming could be written in order to allow interactive 
graphics for design and testing procedures. 
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1 1 — I ( ma ) 




TWO-TERM SINUSOIDAL BASIS 

FIVE-TERM SINUSOIDAL BASIS 

XX X X MACK'S MEASURED DATA 

-M- + 4-+ RESULTS OF THIS STUDY 

FIGURE 7. CURRENT DISTRIBUTION FOR SIMPLE SYMMETRIC 
CENTER-FED DIPOLE 
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MAGNITUDE 


ARGUMENT 


DISTANCE ALONG ELEMENT 


(ma) 


(deg. ) 


FROM CENTER (meters) 


10.069451 76 


-20.91494045546 




.000 


10.31137742 


-23.88716574162 




.025 


10.29630974 


-26.55377940962 




.050 


9.99144960 


-28.92611960589 




.075 


9.37425727 


-30.97158254089 




.100 


8.43725169 


-32.68483783763 




.125 


7.19133077 


-34.07240413432 




.150 


5.66716441 


-35.13865883043 




.175 


3.91449351 


-35.89288284490 




.200 


1 .99942153 


-36.34223309019 




.225 


0.00000000 


-36.49146393058 




.250 


INPUT ADMITTANCE REAL 


IMPUT 


ADMITTANCE IMAGINARY 


(mi 1 isiemens) 






(mil isiemens) 


9.40598988 






-3.59460891 


INPUT IMPEDANCE REAL 


INPUT 


IMPEDANCE IMAGINARY 


(ohms) 






(ohms ) 


92.76685977913 






35.45193909900 


OMEGA RELATIVE VERT. PATTERN 


RELATIVE HORZ. PATTERN 


(deg. ) 


(absolute) 




(absolute) 


000 


1.00000000000 




1 .00000000000 


010 


.97740371048 




1.00000000000 


020 


.91249376631 




1 .00000000000 


030 


.81307695801 




1.00000000000 


040 


.68974122014 




1.00000000000 


050 


.55322994009 




1 .00000000000 


060 


.41223256076 




1 .00000000000 


070 


.27215639036 




1.00000000000 


080 


.13498559880 




1.00000000000 


090 


.00000000000 




1 .00000000000 


100 


.13498559880 




1 .00000000000 


no 


.27215639036 




1 .00000000000 


120 


.41223256076 




1 .00000000000 


130 


.55322994009 




1.00000000000 


140 


.68974122014 




1.00000000000 


150 


.81307695801 




1.00000000000 


160 


.91249376631 




1 .00000000000 


170 


.97740371048 




1.00000000000 


180 


1.00000000000 




1 .00000000000 


The gain of the 


array is 1.64 = 


2. 15db. 



TABLE II. COMPUTER OUTPUT FOR SIMPLE SYMMETRIC 
CENTER-FED DIPOLE 
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APPENDIX A 



GUIDE TO OPERATION OF PROGRAM 

Belov/ is a step by step procedure to follow in order to operate the 
computer program, at the Naval Postgraduate School computer laboratory. 
The software requires the use of the digital computer (XDS9300) and 
the graphics computer (AGT). 

A. FORMATION OF THE DECK 

1. The first card is a 'BOOT' card which is found in the card 
stalls in the laboratory. 

2. The second card is the 'JOB' card which is followed by the 
cards of the program through to the 'DATA' card as shown in the 
enclosed computer program printout. 

3. Following the 'DATA' card is the initial data card which 
contains three digits in the first three columns of the card. The 
first two digits denote the number of arrays to be analyzed (01 to 99). 
The third digit denotes the assigned number of the graphics computer 

to be used (1 or 2) . 

4. After the initial data card, the input data of each of the 
arrays is formed. Each array will have one card plus as many cards as 
there are elements in the array. The first card of each array will 
have the number of elements in the array in the first two columns 

(01 to 10). In the next twenty columns, the value of the wavelength 
in meters is to be entered. 
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5. Following the 'wagelength' card, there is one card for each of 

the elements. The format for each of these cards is as follows. 

Columns 01 to 15: Real value of complex voltage applied to element. 

Example 1.0 

Columns 16 to 30: Imaginary value of complex voltage applied to 

element. Example -1.0 

Columns 31 to 45: Value of the Y-axis, indicating the position of 

the element in meters. Example 0.0 

Columns 46 to 60: Value of half length of element (h) in meters. 

Example 0.25 

Columns 61 to 75: Value of radius of element in meters. 

Example 0.001 

6. As an example, in order to run example one and two of section 



VI of this 


study, the 


cards 


after 


the ‘DATA 


follows : 












Card 1 : 


021 










or 












Card 1: 


022 










Card 2: 


02 


1 .0 








Card 3: 


1 .0 


0.0 


-.005 


0.25 


.000277 


Card 4: 


1.0 


0.0 


.005 


0.25 


.000277 


Card 5: 


01 


1.0 








Card 6 : 


1.0 


0.0 


0.000 


0.25 


.0069 


B. PREPARATION 


OF XDS9300 


DIGITAL 


COMPUTER 



1. Press IDLE switch at the XDS9300 control panel. 

2. Hold the RESET switch depressed and press the POWER switch in 
order to energize the computer. 

3. Turn switch at teletype to the 'on' position. 

4. Press the READY button on the line printer. 
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C. PREPARATION OF THE AGT GRAPHICS COMPUTER 



1. Press ON at the operator's control panel (OCP). 

2. Wait 5 seconds, then press RESET on the OCP. 

3. Press the ON/START switch at the disk drive of the AGT to be 
used. 

4. Press HALT, RESET, and BTSP, in that order on the OCP. 

5. Verify that the 'location counter' (behind the door of the 
left-most rack) indicates zero or all lights out. If not, press 
RESET and BTSP again. 

6. Load a paper tape marked DISK BOOT LOADER (located behind the 
door of the left-most rack) into the tape reader of the AGT's teletype. 

7. Verify that the knob on the front of the AGT’s teletype is in 
the 'line' position. 

8. Put the tape reader switch to 'start.' 

9. Now verify that the 'location counter' indicates 110 (octal). 

If it is not, return to step 1 of this section and begin the 
procedure again. 

10. When step 9 is satisfied, verify that the switches in the 
center of the OCP have the configuration 24000 77776 (octal). 

11. Press HALT, RESET, RUN, and PULSE 1 on the OCP. 

12. If the AGT's teletype does not respond now with a request 

for the operator to input MO/DA/YR, press HALT, RESET, RUN, BUTTON 29, 
and PULSE 1. If teletype does not respond now, return to step 1 of 
this section. 

13. Type RESET( "GATED" ,101 ) ! on the AGT's teletype. 

14. Type GATED! on the AGT's teletype. 

15. Turn vertical and horizontal gain knobs fully clockwise. 
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D. READING IN PROGRAM DECK 

1. Place cards in card hopper. 

2. Push POWER ON button, at the card reader. 

3. Push START button, at the card reader. 

4. Push IDLE, RESET, RUN, and CARDS on the XDS9300 control panel. 
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